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ABSTRACT 

Do spatial distributions of dust grains in galaxies have typical forms, as do spatial distributions 
of stars? We investigate whether or not the distributions resemble uniform foreground screens, as 
commonly assumed by the high-redshift galaxy community. We use rest-frame infrared, ultraviolet, 
and Ha line luminosities of dust-poor and dusty galaxies at z ~ 0 and z ^ 1 to compare measured Ha 
escape fractions with those predicted by the Calzetti attenuation formula. The predictions, based on 
UV escape fractions, overestimate the measured Ha escape fractions for all samples. The interpretation 
of this result for dust-poor z ^ 0 galaxies is that regions with ionizing stars have more dust than 
regions with nonionizing UV-emitting stars. Dust distributions for these galaxies are nonuniform. The 
interpretation of the overestimates for dusty galaxies at both redshifts is less clear. If the attenuation 
formula is inapplicable to these galaxies, perhaps the disagreements are unphysical; perhaps dust 
distributions in these galaxies are uniform. If the attenuation formula does apply, then dusty galaxies 
have nonuniform dust distributions; the distributions are more uniform than they are in dust-poor 
galaxies. A broad range of Ha escape fractions at a given UV escape fraction for z ^ 1 dusty galaxies, 
if real, indicates diverse dust morphologies and the implausibility of the screen assumption. 


1. INTRODUCTION 

A galaxy’s morphology is often defined as its spatial 
distribution of stars. Rest-frame optical images of galax¬ 
ies show diverse arrangements of their stellar compo¬ 
nents: smooth ellipticals, tightly-wound spirals, clumpy 
disks, barred d isks, train-wrecks, fuzzballs . and whatever 
Willman I is (jWillman et al.l I200 !tI 1201 Ill . Little more 
than spatial exte nt is known about distributions of dust 
in du s ty galaxies (lYouneer et al.ll^OlOHDiaz-Santos et al.l 
I20IOI : iSimpson et al.l 1201511 . Yet spatial distributions 
of all galactic components have intrinsic value. If we 
want to learn about galaxies, morphology should be a 
wavelength-independent word. Surprises surely abound: 
for example, galaxies at z 0.3 have dust out to sev¬ 
eral Mpc. Their average dust mass profil e is a constant 
fracti on of th eir average halo mass profile (iMenard et al.l 
I2010L see also iZaritskv 11994) . 

Dust attenuates a galaxy’s intrinsic luminosity at ul¬ 
traviolet (UV) and optical wavelengths. To recover an 
intrinsic luminosity from an emergent luminosity we of¬ 
ten assume that dust grains are distributed as a uni¬ 
formly thick screen between us and the galaxy’s stars 
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(|Calzetti et al.l[r994) . Erom the uniform screen assump¬ 
tion follows an attenuation formula that depends solely 
on wavelength. 

An attenuation formula is denoted fc(A) or k'{X). It 
differs from an extinction formula, which is independent 
of the spatial distribution of dust. When we derive an 
extinction formula, it is valid for the small area surround¬ 
ing a star in the Milky Way, not the large area containing 
a mixture of stars in distant galaxies. 

If we write this equation: 


T / T , _ -|r)-A/2.5 _ 1 Q-E(B-V)fc(A)/2.5 

-^emergent/ -^intrinsic — 

we assume a uniform screen not because k{X) is solely de¬ 
pendent on wavelength, but because ^emergent/^intrinsic 
is solely dependent on wavelength. In an illustrative 
sense ^emergent/^intrinsic Is an attenuation formula con¬ 
volved with a spatial distribution of dust. 

The uniform screen assumption is unrealistic. A low- 
redshift galaxy with a high emergent UV luminosity and 
a low IR luminosity—a dust-poor gala x y— ha s a high 
UV escape fraction. The iCalzetti et al.l I2OOOD formula 
predicts its Ha escape fraction; the prediction overesti¬ 
mates the measured Ha escape fraction. The galaxy’s 
regions with > 10 Mq stars, which ionize gas, have more 
dust than its regions with less massi ve stars that emit in 
the UV and lea ve the gas unionized (jGalzetti et ahlllQ^ 
lCalzettilll997bD . While the discrepancy in amounts of 
dust does not invalidate the screen part of our assump¬ 
tion, it does invalidate the uniform part. 

The screen part of our assumption is suspect, but 
it may be a good approximation to the truth—the al¬ 
ternative assum ption being that dust mixes with stars. 
iLiu et all (I2013D study the dust distribution of M83 at a 
spatial scale of 6 pc. They divide M83 into a center and 
outskirts. Eor each part, they calculate the percentage of 
regions with ionizing stars mixed with dust. Only in the 
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center does this percentage exceed 50. If they degrade 
the spatial resolution to 100 pc, they conclude that the 
dust distribution around ionizing stars is a screen. The 
screen assumption is qualitatively similar to a scenario 
in which dust is dist ributed in a large number of clumps 
(jCalzetti et al.lflo^ . If we choose to think of dust dis¬ 
tributions as clumpy, the number of clumps surrounding 
nonionizing UV-emitting stars in low-redshift dust-poor 
galaxies is 60% of the numb er of clumps surrounding ion¬ 
izing stars (jCalzettilll997all . 

The discrepancy between the uniform assumption and 
real dust distributions in low-redshift dust-poor galaxies 
grows even larger if we consider that galaxies have cen¬ 
ters and outskirts—that galaxies are not point sources. 
A number of studies find that Ha and UV escape frac¬ 
tions increase with increasing radi us (|Boissier et a ll2004 
iPrescott et al.1 l2007t iMunoz-Mateos et al.ll2009l b Ioniz¬ 
ing stars are surrounded by more dust than nonionizing 
UV-emitting stars; ionizing stars at large radii have less 
dust than ionizing stars at small radii, and similarly for 
nonionizing UV-emitting stars. 

We know little about dust distributions in low-redshift 
dusty galaxies, which have high IR luminosities and 
low-to-high emergent UV luminosities. Their IR emit¬ 
ting regions span a range of sizes, from sub- to several 
kpc (iDiaz-Saiitos et al1l2010ll . Measuring a size is hard 
enough; imaging a distribution is even harder. If a dusty 
galaxy’s measured Ha escape fraction agreed with the 
prediction from the lCalzetti et all (|2n0fl[l attenuation for¬ 
mula, we might conclude that the galaxy has a uniform 
dust distribution only if we believe that the prediction 
is valid. If a dusty galaxy’s measured H a escape frac¬ 
tion d isagreed with the prediction from the lCalzetti et all 
(|2000| 1 formula, the galaxy might have a nonuniform dust 
distribution, like that of a dust-poor galaxy. This expla¬ 
nation for the disagreement is not unique. The galaxy 
might have a uniform distribution and obey a different 
fc(A). 

We have nothing more than vague and conflicting 
ideas of dust distribution s in all high-redshift galax¬ 
ies. lOnodera et ^ (j2nifl[ l. iKashino et al.l (j2013h . and 
iPrice et al.l ( 201.^ argue for nonuniform distributions 
in their samples, which c omprise dusty and dust-poor 
galaxies; lErb et akl (|2006ll present evidenc e for uniform 
distrib utions in a similarly mixed sample. iReddv et al.l 
lj2015l l argue that the uniformity of the dust distribution 
depends on the galaxy’s star formation rate: as the star 
formation rate increases the dust distribution becomes 
less uniform. 

In this paper, we ask the question. For dusty galaxies 
at z ^ 0 and 2^1, how do Ha escape fractions r elate to 
the predictions made by the lCalzetti et al.l (|200fl[l atten¬ 
uation formula? We: (1) show that Ha escape fractions 
differ from the prediction; (2) show that the relations 
between Ha and UV escape fractions differ from the re¬ 
lation for low-redshift dust-poor galaxies; and (3) argue 
that an interpretation in the context of dus t distributions 
relies on the shaky assumption that the iCalzetti et al.l 
(j2000f) attenuation formula is universally valid. 

For this paper we assume a cosmology with Hq — 70 
km s“^ Mpc Hin = 0.3, and Ha = 0.7. Our termi¬ 
nology is: (1) a uniform screen is a screen equally thick 
between a region with stars which ionize gas and a region 
with stars which emit in the UV and do not ionize gas; 


(2) a dust-poor galaxy generally has high values of UV 
and Ha escape fractions; and (3) a dusty galaxy is se¬ 
lected from an IR image and generally has low values of 
UV and Ha escape fractions. The samples of dust-poor 
and dusty galaxies are not dichotomous. More details 
are in the following section. 

2. DATA 

Measured and derived quantities are in Tables [2] and 

H 

2.1. Measured quantities 

We use IR, emergent UV, and Ha luminosities to de¬ 
termine escape fractions. In this section we detail our 
catalogs and samples. 

2.1.1. Sample of dusty galaxies at z > 0.7 

Our study uses observations of the GOODS-S, COS¬ 
MOS, and UDS regions. A catalog of Herschel/VKCS 
100 iim sources for GOODS-S comes from iMagnelli et ^ 
dMll; for COSMOS and UDS, we use catalogs pro¬ 
duced for the Herschel survey of CANDELS regions (In- 
ami et ah, in prep). A 100/rm source is defined as a 
> 3cr flux density measurement from PSF fitting to a 
Spitzer/M\FS) 24 /rm source, which is based on an IRAC 
3.6 /rm source prior. In GOODS-S, the catalog has flux 
densities at 24, 70, 100, and 160 /rm; in COSMOS and 
UDS, the catalogs have flux densities at 24, 100, 160, 
250, 350, and 500/rm. All flux densities are > 3cr. We 
do not use the 250, 350, and 500 ^m flux densities in 
the COSMOS and UDS catalogs, to avoid the effects of 
source confusion on PSF fitting. Including them changes 
our results negligibly. 

We associate each 100/xm source with a 1.6/xm (iJ- 
band) source in the c atalogs produced by CANDELS 
(Nay yeri et ah, in prep: lGalametz et al.ll2013l: iGuo et al.1 
[Ml . Each catalog has broadband flux densities at 
many wavelengths, not all of which are common to the 
other two. We require the 1.6/xm source to be a unique 
match within 0.7" of the position of the 100 /xm source’s 
3.6 ^m prior. The match radius is an estimate for max¬ 
imizing the number of unique matches while minimizing 
the number of multiple matches. 

The final quantity of interest for each matched source 
is an Ha line flux or flux limit. These come from our re¬ 
duction of obser vations from the HST /WFC3-IR grism 
survey 3D-HST ([Brammer et al.ll2012[ l. We match each 
100 /xm source to a counterpart in the direct image ac¬ 
companying the grism observations. Two people in our 
group visually inspect the counterpart’s spectrum and 
then review and reconcile any discrepant redshift assign¬ 
ments. The Ha line falls in the grism’s wavelength range 
for galaxies at 0.7 < z < 1.5. The spectra have low 
enough wavelength resolution that the Ha and [N H] 
lines are blended; we assume that the line flux of [Nil] at 
A = 0.6583 /xm is 0.3 times the Ha line flux and that the 
line flux of [N H] at A = 0.6548 /xm is 0.1 times the Ha 
line flux. The ratio of [N H] to Ha line fluxes depends 
on gas-phase metallicity so these values may be incorrect 
for some star-forming galaxies. The values will be correct 
for a significant fraction of galaxies in our sample if the 
sample has a distribution of gas- phase metallicity simila r 
to that of low-redshift galaxies (iKauffmann et alll2003fl . 
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The values will be correct for a larger fraction of galaxies 
in our sample if the sample has a distribution of gas- 
phase metallicity similar to that of low-redshift galaxie s 
which are luminous in the near-IR (| Weiner et al.l[2?)(T^ . 
For each ID spectrum with a visible and uncontaminated 
line complex, we fit a linear model to the continuum and 
Gaussian profiles to the emission lines and extract the Ha 
line flux. The uncertainty on the line flux comes from the 
covariance matrix times the fit’s reduced We replace 
< 5cr Ha line fluxes with limits. We multiply the line 
fluxes and limits by 1.1 as an aperture correction. 

Some sources will be at 0.7 < z < 1.5 but will 
not have visible Ha line flux. To determine line flux 
upper limits for these sources, we match each 100 /rm 
source to a sou rce with a spe c trosco p ic redshift in 


(1200411: IMivnoli et abl (1200, 5n: 

Yamada et al.l ( 

200,511: 

Ravikumar et al. (1200 7t): Iv 

anzella et al.l ( 

200l: 

Lilly et al. (I2009D: iBalestra et al.l (1201011: Fadda et al.l 

(1201011: Simoson et al.l (1201211 

iKurk et al.l (201311: 


et al. (in prep.); and M. Dickinson (private communica¬ 
tion). Operationally both the 100/rm and spectroscopic 
catalogs are matched to the CANDELS catalogs. If 
the 100 /xm source has a spectroscopic redshift that is 
unconfirmed by an uncontaminated grism spectrum, we 
estimate the widths of the invisible Gaussian line profiles 
from the ratio of integrated flux density to peak surface 
brightness for the region of the direct image surrounding 
the galaxy. The direct image is the image taken as part 
of the grism observations. We then fit the same model 
as above to the ID spectrum and extract the Ha line 
flux uncertainty. We use 5 ct limits. We are unable to 
determine limits for galaxies observed by the grism that 
are undetected in the direct image. Most limits are in 
GOODS-S because the spectroscopic completeness is 
higher there than it is in COSMOS or UDS. 

To summarize, our high-redshift sample contains 
galaxies with: (1) detected emission at 100 ^m; (2) de¬ 
tected emission at 1.6/xm; and (3) spectroscopically- 
determined redshifts at 0.7 < z < 1.5. 


2.1.2. Samples of dust-poor and dusty galaxies at z < 0.2 

We use a nonstandard method to determine Ha escape 
fractions, so we test this method first for low-redshift 
dust-poor galaxies which have known dust properties. 
We choose a sample of star-forming galaxies at z ^ 0 
that are analogs of common star-forming ga laxies at high 
redshift fLBAs: lOverzier et al.ll2009L l2011h . 

We use 24/xm flux densities in the Spitzer Herit age 
Archive. Of the 26 sources in lOverzier et al.l (|2011li . 3 
have unreliable 24 /xm flux densities per the recommen¬ 
dations in the enhanced imaging products quick start 
guide; these sources are excluded. The remaining 23 
galaxies have flux densities at 0.15 and 0.23/tm (NUV 
and FUV bands), from GALEX observations. Ha line 
fluxes, from SDSS observations, are from Overzier (pri¬ 
vate communication). We mul tiply the line fluxes by 1.67 
as a fiber-aperture correction (|Overzier et al.l[2(i0^ . 

We also apply our method to a sample of low-redshift 
dusty g alaxies, which have relat ively unknown dust prop¬ 
erties. iHwang fc Gelled (|2013f) collate, for low-redshift 
dust-obscured galaxies (DOGs) and a control sample: 
flux densities at 9, 12, 22, 25, 60, 100, and 140 /tm. 


from IRAS, AKARI, and WISE observations; 0.15 and 
0.23 /xm flux densities, from GALEX observations; and 
Ha line fluxes, from SDSS observations. The fiber- 
aperture corrections for the Ha line fluxes come from 
the differences, in the accompanying 0.62 /xm (r-band) 
images, between Petrosian flux densities and flux densi¬ 
ties in apertures with equal size to the fiber apertures. 

2.2. Derived quantities 

In this section we detail how we estimate total IR and 
emergent UV luminosities, UV continuum power-law in¬ 
dices, and star formation rates for the galaxies in our 
samples. 

2.2.1. IR luminosities 

We estimate a total IR lumin osity (8-1000/xm ; Ltr) 
for each galaxy. We redshift the iCharv &: ElbazI (1200111 
template spectral energy distributions (SEDs) to the dis¬ 
tance of each galaxy. If the galaxy has one measured IR 
flux density, we find the SED that most closely matches 
that flux density and multiply the IR luminosity of the 
SED by the ratio between actual and predicted flux den¬ 
sities. If the galaxy has two or more measured IR flux 
densities, we introduce a multiplicative factor / for each 
SED. We find the best-fit / for each SED; the result 
is an / and value. We choose the SED and its ac¬ 
companying / that results in the minimum of all val¬ 
ues. The galaxy’s Lir is then / times the IR luminosity 
of the SED. In practice, all the high- and low-redshift 
dusty galaxies have two or more measured IR flux densi¬ 
ties, and all the low-redshift dust-poor galaxies have one 
measured IR flux densi ty. 

lOverzi er et al.l (1201111. in their footnote 14, find that 
the iGharv fc ElbazI (l200lh SEDs do not fit low-redshift 
dust-poor galaxies as well as other SEDs. They mea¬ 
sured IR flux density limits at 70 and 160 /xm as well as 
flux densities a t 24 xxm . The footnote explains that the 
IGharv fc Elbad (|20f)lh SEDs, when fit to 24 /xm flux den¬ 
sities, predict flux densities at 70 and 160 /xm that fall 
above their upper limits. Our IR luminosities are high 
compared to theirs—on average, a factor of 1.4 times 
higher; at most, a factor of 2.5 times higher. We present 
results for both sets of IR luminosities for the low-redshift 
dust-poor galaxies. 

IR luminosities are between 5 x 10^° and 1 x IO^^Lq 
for the central 95% of the sample of high-redshift dusty 
galaxies; between 7 x 10^° and 4 x IO^^Lq for the low- 
redshift dusty galaxies; and between 4 x 10^° and 5 x 
IQIILq for the low-redshift dust-poor galaxies. 

2.2.2. UV continuum power-law indices and luminosities 

For galaxies with UV emission from newly formed mas¬ 
sive stars the continuum between 0.125 and 0.263/xm is 
approximately a power law with index /3: 

= CXf^ ( 1 ) 

where S\ is flux density and A is wavelength. We use 
all measured flux densities and uncertainties in this rest- 
frame wavelength range to fit for /3 and the constant 
factor C or for their limits. Measured flux densities may 
be negative—the sky aperture may be brighter than the 
galaxy aperture—which complicates our fitting proce¬ 
dure. When a galaxy has only two measured flux den¬ 
sities in the wavelength range, the solution for /3 and C 
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is analytic, so for a solution to exist the flux densities 
must be positive. If they both are, we make no cuts 
on their signal-to-noise ratio (S/N). If the only positive 
flux density, among multiple measurements, is > 3 (t, we 
calculate a /3 limit using that flux density and a 3cr limit 
for another in the rest-frame wavelength range. Whether 
the limit is an upper or lower one depends on whether 
the detected flux density’s wavelength is lower or higher 
than the limit’s wavelength. We do not estimate a /3 
limit when the only positive flux density, among multi¬ 
ple measurements, is < 3 ct. If more than two measured 
flux densities are positive, regardless of their S/N and 
whether the surplus measurements are negative or posi¬ 
tive, we minimize 

We estimate an emergent UV luminosity or a limit at 
rest-frame 0.16 /rm for each galaxy using its redshift and 
the power-law fit to the rest-frame UV flux densities. We 
refer to the luminosity as 0.16 x Uo.ie- 

2.2.3. Star formation rates 

The equations in iKennicuttI (119981 1 relate IR, emer¬ 
gent UV, and emergent Ha luminosities to star forma¬ 
tion rates (SFRs). We assume that a galaxy’s total SFR; 
(1) is the sum of its IR- and emergent UV-derived SFR; 
and (2) equals the intrinsic (not emergent) Ha-derived 
SFR. The SFR equations are valid under different as¬ 
sumptions for the duration of a galaxy’s star formation 
ep isode. Ou r assum ptions will not conflict with those 
of IKennicuttI (|1998ll if: (1) the SFR has been constant 
for ^ 10® yr, so that the instantaneous Ha-derived SFR 
equals the prolonged UV- and IR-derived SFRs; and (2) 
the dominant source for the IR luminosity is not old stars 
or an active nucleus. 

We must exclude galaxies with active nuclei because 
fractions of their IR, emergent UV, and emergent Ha 
luminosities will be unrelated to newly formed massive 
stars. We identify 5 AGN hosts among t h e high -redshift 
dusty galaxies; we use the iDonlev et al.l (|2012ll criteria, 
which are based on mid-IR flux densities, which all 3 
CANDELS catalogs have. The se criteria will miss some 
AGN hosts (j.Iuneau et al.ll20l'^ . We accept the contam¬ 
inants because in some fields we lack data—for example, 
X-ray data—that we could use to identify the missed 
AGN hosts. 

Among the low-redshift dusty galaxi es, we keep only 
those identihed as star-forming by iHwang fc Gelled 
(j2013f) . They use rest-frame optical line flux ratios to 
discriminate between an AGN host and a star-forming 
galaxy without an AGN. N one of the low-redshi ft dust- 
poor galaxies host an AGN (|Overzier et al.l[2011[) . 

3. RESULTS 

Three star formation rates allow us to compare two at¬ 
tenuation values. We eschew magnitude attenuation in 
favor of the escape fraction, which is the ratio of emergent 
to intrinsic luminosity. We define the following quanti¬ 
ties. 

• fesc(0.16p,m) is the escape fraction at rest-frame 
0.16^m. In relation to AyV) fesc(0.16/im) = 

10“Auv/2.5 

• fesc(Ha) is the escape fraction of the Ha line lumi¬ 
nosity. 


• fesc(0.66/rm) is the escape fraction extrapo¬ 
lated from fesc(0.16/xm) to the wavelength of 
Ha—the attenuation form ula prediction . Fo r 
the attenuation formula in iGalzetti et al.l (j200(lfl , 

fesc(0.66/xm) = fesc(0.16/xm)°®®. 

IGalzetti et al.l (|1994ll argue that, in low-redshift dust- 
poor galaxies, regions with ionizing stars are spatially 
separate from regions with nonionizing UV-emitting 
stars. In this framework fesc(0.66 /xm) is the escape frac¬ 
tion of the 0.66 /im continuum emission from low-to-high 
mass stars and fesc(HQ:) is the escape fraction of the 
0.66 /im line emission from gas ionized by high mass stars. 
These galaxies have fesc(Ha) = fesc (0-16/im)° ®^, equiv¬ 
alent to UsnfHg ) = fesc(0.66/im)^-^ (|Galzetti et al.lHO^ 
lCalzettilll997bf) . Because the ionizing and nonionizing 
stars are spatially separate, the escape fractions disagree 
because the dust distributions are nonuniform. The ion¬ 
izing stars are surrounded by more dust than are the 
nonionizing stars. 

The Ha escape fractions in IGalzetti et al.l (119941) 
come from the ratios of Ha to H,5 line luminosi¬ 
ties. In the absence of dust the ratio is 2.86, so 
(THa/TH/3)/2.86 = fesc(Ha)/fesc(H/3). Since fesc(H;9) is 
a power-law function of fesc(IIa), the ratio of line lumi¬ 
nosities uniquely determines fesc(IIa). 

The UV and 0.66/im escape fractions in iCalzetti et al.l 
(| 199411 come from comparing spectra of galaxies with Ha 
escape fractions of 1 to spectra of galaxies with Ha es¬ 
cape fractions less than 1. The Ha and H/3 lines in the 
spectra must be masked or replaced with continuum val¬ 
ues. The literat ure relat i on is t echnically fesc(IIa) = 
fesc(0.66/im)^ ®. iGalzettil (I1997U1 uses the Milky Way 
extinction formula of lHowarthI ( 19831 1 to determine Ha 
escape fractions for low-redshift dust-poor galaxies; if we 
instead use the attenuation formula, we get fesc(IIa) = 

fesc(0.66/im)2U 

Our hypothesis is that the Ha and 0.66 /xm escape 
fractions agree—that dusty galaxies have uniform dust 
distributions. Our method for determining Ha escape 
fractions is nonstandard; we lack H/3 luminosities for the 
z ^ I dusty galaxies. We first apply it to a sample of low- 
red shift dus t -poor g alaxies with the aim of reproducing 
the IGalzettil (jl997bl l result. 

3.1. Reproducing the relation between escape fractions 
for dust-poor galaxies 

The sample, from lOverzier et al.l (|2011ll . contains 
galaxies that lie roughly on the low-redshift IRX-/3 re¬ 
lation. This relation is a blunt instrument we use to 
determine whether or not galaxies obey an attenuation 
formula at UV wavelengths; fesc(0.16/xm) is a function of 
IRX and an attenuation formula relates fesc(0.16/xm) to 
p. Since these galaxies lie on the relation, we expect them 
to have dust propertie s that conform with t he properties 
of the galaxies in the IGalzetti et al.l (|200nr ) sample. We 
calculate fesc(IIa) and fesc(0.16/xm) for each galaxy in 
the following way. 

• fesc(0.16/xm) = (Lir/(1.68 x 0.16 x Lo,i6)+l)“^- 

This is equivalent to fesc(0.16/xm) = (emergent UV 
SFR)/(IR SFR -L emergent UV SFR). 

• Following the assumptions in 12.2.31 fesc(IIa) = 
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(emergent Ha SFR)/(IR SFR + emergent UV 
SFR). 

We use these quantities to find, for the sample, a best-fit 
power-law index q for the form fesc(Ha) = fesc(0-16 ^m)'^. 

The top left panel of Fig. [T] shows that fesc(Ha) = 
fesc(0.16for the low-redshift dust -poor galaxies. 
If we use the IR and UV luminosities from lOverzier et "m 
I 2 OIII fesc(Ha) = fesc(0.16/rm)°'®'*. We lack uncertain¬ 
ties on any SFRs so we also perform a robust linear re¬ 
gression, which requires a fit to the form logfesc(Ha) = 
q logfesc(0.16 ^m). It returns a similar result, as does a 
bootstrap resam pling. We recover a relation close to the 
ICalzettil (|1997bf) relation for similar galaxies, so we have 
some confidence that our procedure is grossly correct. 
We stress, again, that there are no evolutionary connec¬ 
tions between the low- and h igh-redshift gala xy samples; 
we aim only to reproduce the lCalzettil (Il997bll result with 
our nonstandard method of determining fesc(Ha). 

3.2. A new relation for dusty galaxies at low redshift? 

The top right panel of Fig. [T] shows that fesc(Ha) = 
fesc(0.16for the low-redshift dusty galaxies. We 
use uncertainties only on the emergent Ha-derived 
SFR—we treat it as the dependent variable—in the 
weighted fit. The relation for dust-poor galaxies is seem¬ 
ingly inapplicable to dusty galaxies at low redshift. In 
il3.4l andHlwe explain why we preface ‘inapplicable’ with 
‘seemingly.’ 

3.3. A corresponding relation for dusty galaxies at high 

redshift? 

The bottom left panel of Fig. [T]shows that fesc(Ha) = 
fesc(0.16^m)°-®® for high-redshift dusty galaxies. We 
leave limits out of this fit. The uncertainty in the in¬ 
dex is very small; this is not surprising given the number 
of galaxies. The dispersion about the best-fit relation is 
large; the central 50% of galaxies with Ha-detected SFRs 
lies between fesc(Ha) = fesc(0.16/im)°-^® and fesc(Ha) = 
fesc(0.16A survival regression, which incorpo¬ 
rates Ha SFR limits in a fit to the form logfesc(Ha) = 
q logfesc(0.16^m), finds q = 0.47, as does a simple linear 
regression using only the Ha-detected SFRs. None of the 
galaxies have limits for both the UV and Ha SFRs. The 
main effect of including limits in the regression is to in¬ 
crease the standard error of the sample’s best-fit index; 
a marginally different index is a secondary effect. Table 
[T] summarizes the relevant quantities for the main fit to 
each sample of galaxies. 

3.4. Issues with the regressions 

The escape fraction relations of dusty galaxies are sub¬ 
ject to selection effects. The best-fit relation for the low- 
redshift dusty galaxies seems to: (1) underestimate Ha 
escape fractions at low UV escape fractions; and (2) over¬ 
estimate Ha escape fractions at high UV escape fractions 
(top right panel of Fig. [T]). This may reflect our con¬ 
straint that fesc(Ha) =1 when fesc(0.16 /im) =1, but the 
lack of galaxies with Ha escape fractions at or below 
the relation may be due to galaxies with Ha SFRs be¬ 
low the limit of the sample. We do not measure Ha line 
fl ux limits for galaxies th at meet the selection criteria 
of lHwang fc Gell^ l|2013ll and have undetected Ha line 


fluxes. If we replace all Ha-detected SFRs with SFRs 
corresponding to the minimum Ha line flux of the sam¬ 
ple, we fill in the space around the relation. 

High-redshift dusty galaxies with UV escape fractions 
< 10“^ seem to have less dispersion in their Ha es¬ 
cape fractions than do galaxies with UV escape fractions 
> 10“^. Despite our best efforts to measure Ha line flux 
limits, we still lack them for galaxies lacking spectro¬ 
scopic redshifts. If we replace all Ha-detected SFRs with 
SFRs corresponding to a reasonable Ha line flux limit for 
the grism observations (5 x 10“^^ erg s“^ cm“^), we fill 
in the space with low Ha escape fractions at low UV es¬ 
cape fractions. The decrease in dispersion is probably 
an artifact of the grism detection limit and a paucity of 
spectroscopic redshifts in UDS and the part of COSMOS 
in CANDELS. 


4. DISCUSSION 


If we assume th at dusty galaxies: (1) obey the 
iCalzetti et al.l (|200n[ l attenuation formula; and (2) were 
to all have the same UV spectrum in the absence of 
dust; our results indicate that their Ha and 0.66/rm es¬ 
cape fractions disagree. Unlike low-redshift dust-poor 
galaxies, dusty galaxies have fesc(Ha) = fesc(0.66/rm)^-®; 
like low-redshift dust-poor galaxies, dusty galaxies have 
nonuniform dust distributions. The relations between 
fBar.(Ha) and L s n(0-66 um) are similar to th o se found by 
loirodera et al.f (1201011.iKashino et all (l2013ll .[P rice et all 
( 201311. and" Reddv et al.l ( 201511 and conflict with those 
of lErb et al.l (1200611. These studie s assume the applica¬ 
bility of the ICalzetti et al.l ()2000l l attenuation formula, 
and use samples of high-redshift galaxies that are not 
necessarily dusty. Their UV escape fractions come from 
stellar population modeling, while ours come from IRX; 
their Ha escape fractions come from ratios of averaged 
Ha to averaged H/3 luminosities, or ratios of Ha to UV 
luminosities, while ours come from assumptions about 
SFRs. 

However, d usty galaxies may disobey the 
ICalzetti et al.l (1200011 attenuation formula, which 
complicates our interpretation of the more basic rela¬ 
tion between UV and 0.66 /rm escape fractions—the 


e mitting stars JII; 

ri [Miif[2(M E 


Wvi IRuat et all __ ______ 

Kriek &: ConrovI 1201311 . The true attenuation formula 


(IReddv et al.j 

201 c 

IPenner et al.l 

2012 


for dusty galaxies may have a higher 0.66/rm escape 
fraction for a given UV escape fraction than does the 
ICalzetti et al.l (1200011 formula. A formula modified ac- 
cording to iKriek fc Conro^ l|2013l l with 6 = —0.2 and no 
bump at 0.2175 ^m has fesc(0.66/xm) = fesc(0.16 
instead of fssr(0- 66/xm) = fBsr,(0-1 6 The com¬ 

bination of this IKriek &: ConrovI (|2013ll formula and 


our r esult implies that fesc(Ha) = fesc(0.66/xm)^'^. 
The (|Reddv et al.l 1201511 attenuation formula has 

fesc(0.66/xm) = fesc(0.16/xm)°-^3, implying that 

fesc(Ha) = fesc(0.66/xm)^'^. If the true formula in¬ 
stead has a much lower 0.66 /xm escape fraction for a 
given UV escape fraction than does the ICalzetti et al.1 
(|2000ll formula—say, fesc(0.66^m) = fesc(0.16/xm)°-^® — 
we conclude that for dusty galaxies the Ha and 0.66 ^m 
escape fractions agree. Our conclusion regarding the 
uniformity of dust distributions depends on the assumed 
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Table 1 

Summary of regressions 


Sample 

<? 

Std. error of q 

Residual std. error 

DOF 

Low-z dust-poor 

0.799 

0.094 

6.0 

13 

Low-z dusty 

0.548 

0.029 

67.0 

54 

High-z dusty 

0.549 

0.008 

6.6 

320 


g Au^(mag) ^ 



O) 

03 

E 


f,3,(0.16|im) (L. 


■emergen/intrinsic) 


Figure 1. Top left: The Ha escape fraction, fesc(HQ), as a function of the UV escape fraction, feac(0.16 /rm), for the low- 
redshift dust-po or galaxies. The das hed purple line shows fesc(0.66 /im), the extrapolation of fesc(0.16 /rm) to the wavelength 
of Ha using the lCalzetti et ahl Il2000f) att enuation f ormula , as a function of fesc(0.16 fim). This is the uniform dust distribution 
hypothesis for this attenuation formula. ICalzettil (|1997bf l. who study another sample of low-redshift dust-poor galaxies, find 
that Ha escape fractions are lower than this extrapolation (dotted black line). We find a similar result (solid orange line). We 
have some confidence that our Ha escape fractions are grossly correct. Top right: fesc(Ha) as a function of fesc(0.16 /rm) for 
the low-redshift dusty galaxies. In low-redshift dusty galaxies, regions with ionizing stars have lower escape fractions than the 
extrapolation and higher escape fractions than they would have in low-redshift dust-poor galaxies. Bottom left: fesc(Ha) as 
a function of fesc(0.16 /rm) for the high-redshift dusty galaxies. In high-redshift dusty galaxies, regions with ionizing stars 
have lower escape fractions than the extrapolation and higher escape fractions than they would have in low-redshift dust-poor 
galaxies. A conclusion about the uniformity of dust distributions in high-redshift dusty galaxies is subject to an assumption 
about the attenuation formula. 

atte nuation formul a . We cannot reject our hypothesis. galaxies (bottom left panel of Fig. [Ij, so if it is the 

The iCalzetti et al.l ()2000( ) formula does demarcate the true formula then at least some dusty galaxies have 

upper boundary of fesc(Hct) for the high-redshift dusty uniform dust distributions. The unambiguous statement 
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we make is that dusty galaxies have relations between 
Ha and UV escape fractions that are different from the 
relation for low-redshift dust-poor galaxies. 

The dispersion about the best-fit relation for high- 
redshift dusty galaxies is large, yet it is susceptible to 
our possibly invalid assumption that a galaxy’s instanta¬ 
neous SFR has, for at least 10® yr, equaled its prolonged 
SFR. For example, if most massive stars formed over 15 
Myr instead of 300 Myr, we should multiply To.ie and 
Lir by larger numbers than we do here to derive the UV 
and IR SFRs. The true UV SFR is 57% higher and th e 
true IR SFR is 63% higher (|Madau fc pickinsonll2014D . 
The true UV escape fraction is little different; because 
the true total SFR is higher, the true Ha escape fraction 
is lower. Another scenario is that massive stars stopped 
forming within the last 15 Myr. Perhaps some of the 
high-redshift dusty galaxies have complicated star for¬ 
mation histories; the combined effect of deriving correct 
SFRs might lead to lower dispersion. 

There are few ways to resolve the uncertainties re¬ 
garding st ar formation h i storie s and the applicabil¬ 
ity of the ICalzetti et al.l (|2000f) attenuation formula 
to dusty galaxies. One (partial) way is to use the 
formula and fesc(Ha) to predict the escape fraction 
of the H/3 luminosity, fesc(H/3). The ratio of in¬ 
trinsic Ha to H/3 luminosities is 2.86 under reason¬ 
able assumptions; with H/3 observations, we can use 
(LHa/iH/3)/2.86 = fesc(Ha)/fesc(H/3) and from this de¬ 
termine fesc(H/3). If the predicted and determined H/3 es¬ 
cape fractions agree, the assumed attenuation formula is 
probably valid at optical wavelengths; the assumed star 
formation histories are also probably valid. The histo¬ 
ries are valid because the prediction is history-dependent 
and the determination is history-independent—both line 
luminosities are from ionized gas surrounding the same 
stars. If the predicted and determined fractions disagree, 
the attenuation formula at optical wavelengths can be 
different for dusty galaxies or their true fesc(Ha) can be 
different. The formula can differ due to dust not dis¬ 
tributed a s a uniform fo reground screen; see table 1 and 
fig. 3b of iCalzettil 120011 Fig. [2] shows this test applied 
to the low-redshift dusty galaxies. To apply this test 
to high-redshift dusty galaxies, we need near-IR spectra 
that cover the H/3 line. If the dispersion of Ha escape 
fractions at a given UV escape fraction for high-redshift 
dusty galaxies remains, there must be a diversity of dust 
distributions. Either way, this test provides no infor¬ 
mation on the relation between UV and 0.66 /im escape 
fractions; more extrapolations are necessary to defini¬ 
tively interpret Fig. [T] 

5. SUMMARY 

1. We measure Ha line luminosities and limits, from 
HST /WFC3-IR grism spectra, for z ^ 1 galaxies 
with detected 100 fim emission. We determine rest- 
frame UV continuum power-law indices {f3 values) 
and ratios of IR to UV luminosities (IRX values). 

2. For each galaxy we determine fesc(0.16 /rm), the es¬ 
cape fraction at 0.16 /im, from IRX. We determine 
fesc(Ha), the Ha escape fraction, from the Ha lu¬ 
minosity and the total star formation rate. 

3. For the lOverzier et al.l l|2011f) sample of low- 
redshift dust-poor galaxies, we recover the ICalzettil 



^esc (Hot) (Lgn^ergen/^ intrinsic) 


Figure 2. H/3 escape fraction, fesc(H/3), as a function of Ho; es¬ 
cape fraction, fesc(HQ:), for the low-redshift dusty galaxies. We 
determine fesc(H;d ) using fesc(HQ ) and the ratios of Ha and H/3 
luminosities. The ICalzetti et al.l l|200Qll attenuation formula has 
fesc(H/3) = fesc(Ha)^'^® (dashed purple line). If H/3 measurements 
of high-redshift dusty galaxies agree with the prediction, we can 
resolve uncer tainties related to sta r formation histories and the va¬ 
lidity of the ICalzetti et al.l I|2Q0QI1 attenuation formula at optical 
wavelengths; otherwise, we might not. 

(|1997bf) relation between Ha and UV escape frac¬ 
tions: fesc(Ha) = fesc(0.16/tm)° ®°. Our nonstan¬ 
dard method for determining Ha escape fractions 
should be grossly correct. 

4. For the iHwang fc Gelleil (|2013[ 1 sample of low- 
redshift dusty galaxies, we find that fesc(Ha) = 

fesc( 0 . 16 /rm) 0 ®®. 

5. For our sample of dusty galaxies at 0.7 < z < 1.5, 
we find that fesc(Ha) = fesc(0.16/tm)° ®®. The dis¬ 
persion about the best-fit relation is large. 

6. The interpretation of the results fo r dusty galax¬ 
i es is unclear. If we assume that the ICalzetti et ahl 
(j2000l l attenuation formula applies to these galax¬ 
i es, ou r results agree with th o se of lOnodera et al.l 
(1201011. iKashino et al.l (2(H^ . iPrice et al.l (j2013ir 
an d iReddv et al.l (12015 1 and conflict with those 
of lErb et al.l (|2006f l: fesc(Ha) = fesc(0.66/im)^ ®. 
They study samples comprising dust-poor and 
dusty galaxies. Dusty galaxies have nonuniform 
dust distrib utions. However, du sty galaxies may 
disobey the ICalzetti et al.l (1200011 attenuation for¬ 
mula. 

7. Measurements of the ratios of Ha to H/3 luminosi¬ 
ties may decrease the dispersion in the relation be¬ 
tween Ha and UV escape fractions for dusty galax¬ 
ies at z ~ 1. If the dispersion is real, these galaxies 
have diverse dust distributions. 

ALMA is beginning to produce images of the dust dis¬ 
tributions in high-redshift galaxies. These images lead to 






















































direct tests, independent of the attenuation formula, of 
the uniform screen assumption. They also have the po¬ 
tential to reveal the parts of galaxies that may never be 
detected at UV and optical wavelengths: surprises surely 
abound. 
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Table 2 

Measured and derived quantities for the sample at z > 0.7 


RA, IR 

Dec, IR 

RA, opt. 

Dec, opt. 

z 


era 0.16 X Lo.16 


Lua 


Lho lim. 

Lm, 


deg 

deg 

deg 

deg 



Lq 

L© 

Lq 

L© 

L© 

Lq 

Lq 


GOODS-S 


53.113487 

-27.933195 

53.113470 

-27.933294 

1.098° 

0.53 

0.34 

l.le-l-10 

4.9e-|-08 

5.6e-|-08 

2.1e-|-07 .. . 

6.2e-|-ll 

1.9e-|-10 

53.096371 

-27.925898 

53.096385 

-27.925972 

1.146 

-1.89 

0.41 

7.9e-|-09 

3-3e-|-08 

3.3e-|-08 

2.6e-|-07 .. . 

2.2e-|-ll 

4.9e-|-10 

53.090958 

-27.922218 

53.090977 

-27.922288 

1.390 

1.41 

0.48 

8.4e-|-08 

4.5e-|-08 

3.1C-I-08 

2.2e-|-07 ... 

5.5e-|-ll 

2.7e-|-10 

53.184750 

-27.920420 

53.184801 

-27.920432 

0.953 

-1.68 

0.15 

1.2e-|-10 

3-9e-|-08 

1.6e-|-08 

2.3e-|-07 ... 

9.3e-|-10 

1.7e-|-07 

53.104836 

-27.913826 

53.104830 

-27.913926 

1.090 

0.61 

0.61 

3.5e-|-09 

2.7e-|-08 

1.2e-|-08 

1.7e-|-07 ... 

1.8e-|-ll 

6.2e-|-09 


Note. — Columns: RA and Dec, IR, are positions in the Herschel catalogs. They are based on positions of Spitzer/IUAC priors. RA and Dec, opt., are positions in the CANDELS catalogs. 
Lhq lim. is a 5cr upper limit to the Ho line luminosity. The entire table is published in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 

Q 

Spectroscopic redshift solely from grism observations. 






Table 3 

Measured and derived quantities for the samples at 2 < 0.2 


RA 

Dec 

z 

/3 

<7/3 

0.16 X 1/0.16 



Lua <7Lho 

^IR 


deg 

deg 




Lq 

L© 

L© Lq 

L© Lq 

Lq 

Lq 


Dust-poor galaxies 

30.987125 

-8.132919 

0.189 

-1.66 

0.04 

3.5e-|-10 

2.9e+08 

1.4e-|-08 ... 

4.6e+08 ... 

8.5e-|-10 

1.8e-|-09 

33.452250 

12.997628 

0.219 

-0.45 

0.05 

4.3e-|-10 

4.3e+08 

3.6e-|-07 ... 

1.4e+08 ... 

5.4e-|-ll 

3-0e-l-09 

52.191625 

1.197458 

0.142 

-1.68 

0.04 

1.9e-|-10 

1.7e+08 

7.4e-|-07 ... 

2.4e+08 ... 

5.4e-|-10 

l.Oe-l-09 

59.391667 

-5.622139 

0.204 

-1.28 

0.06 

2.9e-|-10 

3.5e+08 

9.4e-|-07 ... 

3.5e+08 ... 

l.le-l-11 

1.9e-|-09 

125.007167 

50.844211 

0.217 

-1.55 

0.08 

2.6e-|-10 

4.0e+08 

2.4e-|-08 ... 

8.1e+08 ... 

4.8e-|-ll 

4.9e-|-09 


Note. 


The entire table is published in the electronic edition of ApJ. A portion is shown here for guidance regarding its form and content. 









